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ABSTRACT 

We present a simple, semianalytical approach to compute the line emission 
of C"*", C° and CO in Photodissociation Regions of varying metallicity that 
depends on very few parameters and naturally incorporates the dumpiness of 
the interstellar medium into the calculations. We use it to study the dependence 
of the I[cii]/Ico and I[ci]/Ico hne ratios on the metallicity. We show how to 
modify it to include the effects of density and radiation field, and we have 
compared it with observational data. We find that the model explains the 
observed trend of enhanced [C II] /CO line ratio with decreasing metallicity 
as the natural result of the augmented fraction of photodissociated gas in a 
clump. We show that enhanced [C II]/CO ratios can be produced by lowering 
the upper limit of the clump size distribution, as may happen in 30 Doradus. 
We also find that the available data favors a [C I] /CO intensity ratio essentially 
independent of metallicity, albeit the paucity of observations does not exclude 
other possibilities. This is difficult to understand if most of the C° is produced 
by UV photons as is the case with C+. Finally, we study the prediction of 
the model for the trend of the Xqo factor with metallicity. Comparison with 
previous observational studies yields good agreement. 



Subject headings: ISM: abundances — ISM: general — galaxies: ISM — radio 
lines: ISM 



1. Introduction 

1.1. Photodissociation regions and the metallicity of the ISM 

A photodissociation region (PDR or photon dominated region) is that portion of the 
interstellar medium where ultraviolet (UV) radiation dominates the physical and chemical 
processes. The first theoretical studies of PDRs (e.g., Tielens and Hollenbach 1985a; 
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Wolfire, Tielens, & HoUenbach 1990) considered the strong UV fields in the vicinity of OB 
stars. These models and others have shown that UV photons not energetic enough to ionize 
hydrogen [hiy > 13.6 eV) escape the H II region and interact with the surrounding molecular 
gas. These soft UV photons photodissociate molecules and ionize atoms with ionization 
potentials less than 13.6 eV, such as carbon, silicon, and sulfur, creating a transition zone 
between the ionized and molecular gas. 

Theoretical PDR models (e.g., Tielens and Hollenbach 1985a) predict that the primary 
coolants of PDRs arc far-infrared (FIR) fine-structure atomic lines, particularly those of 
[O 1] and [C 11]. Hence, these lines arc expected to be quite bright near regions of massive 
star formation. This expectation is borne out by observation: the [C II] 158 fim line alone 
accounts for 0.1 to 1% of the total far-infrared luminosity in spiral galaxies (Stacey et al. 
1991), corresponding to a luminosity of up to ~ 10^ Lq. The [C II] line is by far the most 
luminous spectral fine in galaxies. Because the FIR fines are so bright, they are often used 
to infer the physical conditions in gas near recently formed massive stars. 

Wolfire et al. (1990) suggest that the most important observational diagnostics for 
PDRs are: (1) CO, which provides an estimate of the molecular gas mass, (2) [C II] and 
[O I] FIR fine structure line emission, which provide information on the densities and 
temperatures, and (3) FIR continuum emission, which, in combination with (1) and (2) 
provide the incident UV flux, the stellar luminosities, and the dust opacities. Recently, 
[C I] submillimeter line emission has also been shown to be an important diagnostic of 
the PDR/molecular cloud interface (e.g.. Plume, Jaffc, & Kccnc 1994; Kccnc et al. 1993), 
especially in translucent clouds with Ay ~ 1 — 2 (Ingalls, Bania, & Jackson 1994 ; Ingalls 
et al. 1997). 

Although there has been considerable progress in the study of PDRs, it is still difficult 
to separate the influences of various physical parameters such as UV fleld strength and 
metallicity on the structure and composition of the interstellar medium. In this paper 
we examine the effects of metalficity on the structure of PDRs and the intensity of their 
diagnostic lines. 

Both theory and observation suggest that metalficity profoundly alters the properties 
and structure of the photodissociation regions (e.g. Israel 1988; Maloney 1990; Madden et 
al. 1997). Several factors contribute to the differences between high and low metalficity 
systems. Because low metallicity systems have fewer heavy elements, the gas phase 
chemistry is altered, resulting in different equilibrium abundances of the molecular species 
(e.g., van Dishoeck & Black 1988; Lequeux et al. 1994). Grain surface chemistry is also 
altered because a deficiency of heavy elements (especially Si, C and S) lowers the dust-to-gas 
ratio. Finally, since there is less dust shielding, UV photons penetrate more deeply into 
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molecular clouds. Molecular hydrogen is largely unaffected by this enhanced UV radiation 
because of the strong H2 self- shielding in the Werner and Lyman bands and the mutual 
shielding between coincident lines of H and H2 (Abgrall ct al. 1992). However, molecules 
which are not so strongly self-shielding, such as CO, are destroyed via photodissocation 
everywhere except in the most opaque clumps (Maloney & Black 1988, Elmegreen 1989). 
Thus, a low- met alii city system contains large regions where hydrogen remains molecular 
but the usual tracers of molecular gas hke CO are absent. Accordingly, molecular gas in 
low-metaUicity systems is not well traced by CO. 

The paucity of CO emission from molecular gas in low-metallicity systems is confirmed 
by observations of the ISM in metal-poor dwarf galaxies. Dwarf galaxies systematically 
display very low levels of CO emission (e.g., Tacconi & Young 1987). CO observations of 
the Large Magellanic Cloud (LMC, Cohen et al. 1988; Israel et al. 1993) and the Small 
Magellanic Cloud (SMC, Rubio, Lequeux, & Boulanger 1993) show that the ratio between 
the CO luminosity and the virial mass is much smaller in LMC and SMC clouds than that in 
Galactic clouds. Furthermore, a comparison of outer Galaxy clouds with those in the SMC 
(where the metallicity is 1/6 relative to the Orion Nebula) shows that CO is less luminous 
by a factor of 10 — 20 over large spatial scales (~ 200 pc) (Rubio et al. 1993). Wilson (1995) 
and Arimoto, Sbfue, & Tsujimoto (1996) have examined the value of the Xco factor, the 
ratio between the CO luminosity and the hydrogen column density, in galaxies spanning 
a large range of metaUicities. They find a strong dependence of Xqo on metallicity: for a 
given molecular gas mass, lower metaUicity systems emit less CO J = 1 — > 0. These results 
all suggest that CO is fainter in low-metallicity clouds due to enhanced photodissociation. 

Observations of the [C II] 158 jim line also provide strong evidence that photodissocation 
regions are greatly enhanced in dwarf galaxies. Ionized carbon is the end product of the 
photodissociation of CO. Recent KAO observations of the [C II] 158 iim transition in the 
LMC (Poglitsch et al. 1995, Israel et al. 1996) and in IC 10 (Madden et al. 1997) show 
that the [C II] emission is extended and extremely bright. In fact, the ratio of [C II]/FIR 
continuum emission in the LMC is about 1%, much higher than in the Galaxy and most 
galactic nuclei. Furthermore, the [C 11] /CO luminosity ratios for H 11 regions in the LMC 
arc larger than those in the Milky Way by a factor of ~ 5 for N160 and a factor of ~ 50 
for 30 Doradus (Israel et al. 1996). Similar results are found for the northern dwarf 
galaxy IC 10 (Madden et al. 1997). The dense ISM in dwarf galaxies seems to be a vast 
photodissociation region, where CO is present only in well shielded clumps. 

To better understand the transition between ionized, neutral atomic, and molecular 
forms of carbon in low-metallicity dwarf galaxies, we have begun a project to study their 
submillimeter [C I] lines. Neutral carbon emission was detected from a low-metallicity 
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system for the first time using the Antarctic SubmiUimeter Telescope and Remote 
Observatory (AST/RO), a new 1.7 m diameter telescope recently deployed at the SoTith 
Pole (Stark et al. 1998). The [C I] line has been measured measured towards two positions 
in the LMC (Stark et al. 1997): N 159, the star forming region near one of the brightest 
CO J = 1 ^ peaks in the LMC, and 30 Doradus, the brightest H 11 region in the local 
Universe . We have also detected [C 1] emission from the dwarf galaxy IC 10 using the 
JCMT at Hawaii (Bolatto et al. 1998). 

These new data prompt us to model the intensities of the various forms of carbon in 
molecular clouds of varying metallicity. 

1.2. Modelling low metallicity PDRs 

There are several models of photodissociation regions that treat the chemistry, radiative 
transfer and energy balance of molecular clouds in a detailed, self consistent and highly 
sophisticated way (e.g.: Tielcns & Hollenbach 1985a; van Dishoeck & Black 1988; Le 
Bourlot et al. 1993; Sternberg & Dalgarno 1995; Spaans 1996; Pak et al. 1998). Many 
of them perform the calculations in a plane-parallel, uniform density geometry. However, 
molecular clouds are clumpy objects, highly permeable to the UV radiation that drives 
most of the chemistry in the PDR (e.g., Stutzki et al. 1988). Applying the results of 
plane-parallel calculations in the context of clumpy PDRs (e.g., Poglitsch et al. 1995; 
Madden et al. 1997) necessitates the additional step of assuming different beam filling 
fractions for different species; otherwise the observed line ratios of the different carbon 
species cannot be explained. This problem is even more acute for low metallicity systems, 
where the UV field permeates the clouds even more pervasively, due to the higher gas to 
dust ratios and the smaller dust shielding. Spaans and coworkers (e.g., Spaans 1996; Spaans 
& van Dishoeck 1997; Spaans & Neufeld 1997) have made numerical 2-D and 3-D Monte 
Carlo models of PDRs that incorporate the effects of dumpiness in a very complete way, 
albeit only for particular realizations of the distribution of clumps within the cloud. 

The purpose of the model presented here is to bridge the gap between the plane-parallel, 
uniform slab calculations and the numerical models of clumpy PDRs. This model naturally 
incorporates the effects of different beam fiUing fractions for the different carbon species, 
and the predicted line ratios depend only on a few free parameters that can be directly 
linked to observations. We do not intend to compute accurately the chemistry, radiative 
transfer and energy balance. Instead, we use simple fits to the chemistry calculations of 
previous plane-parallel models to include the effects of density and radiation fields in our 
calculations, and we rely on observations to determine the values of our model parameters. 
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In §0 we will state some of our initial simplifying assumptions, in §^ we will develop the 
basic mathematics of the model and compute the column densities of the carbon species, in 
^ we will compute the line intensity ratios, in §^ we will discuss the values for the model 
parameters, in §^ we will discuss how these values are affected by changes in the gas density 
and radiation field, in §0 we will compare our model results with the available observations, 
and in §H we will discuss the relationship between this model and the observed CO intensity 
to H2 column density conversion factor. 

2. The model 

We assume that the gas phase of the ISM is constituted by uniform density, unresolved, 
spherical clumps of radius R immersed in an isotropic UV field. The assumption of an 
isotropic UV field is justified in view of the large UV back-scattering coefficient for dust 
grains (e.g.. Savage & Mathis 1979). We consider a distribution of clump radii in the 
telescope beam of the form dN/dR = R~^. This type of power law for the distribution 
of clump sizes has been observed in molecular clouds with exponent /? ~ 3.3 (Stutzki & 
Glisten 1990; Elmegreen & Falgarone 1996). To simplify the mathematics of the model 
we will take this exponent to be /3 = 3. Changes of order 10% in f3 affect minimally the 
model results. The model implicitly considers unresolved clumps; our resolution element 
is greater than the maximum clump size. The observed intensities then result from the 
combined contributions of the ensemble of clumps within the telescope beam. We use the 
same average density n for all the clumps. Studies of the mass and size distributions of 
clumps in Milky Way molecular clouds (Leisawitz 1990; Elmegreen & Falgarone 1996 and 
references therein) reveal that the clump masses follow a power law M oc i.e., the 

clump volume density n appears to be roughly constant. Thus, we are not incorporating a 
density distribution or a size-density relationship (Larson 1981; Myers 1983). Throughout 
this work we assume that the dust-to-gas ratio is proportional to the metallicity of the 
interstellar medium. 

According to PDR theory, a typical clump has three well defined regions: an outermost 
zone where C^ is the dominant form of carbon, a core where almost all the carbon is locked 
in CO molecules, and a transition region in between where most of the carbon is in neutral 
atomic form (C). We model each clump as consisting of three concentric zones with sharp 
boundaries, where all the carbon is in the dominant form of the species corresponding to 
each zone. 
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3. The column density of the carbon species 



3.1. Model A 



Let the physical extent of the C*^ region be h. We assume that the following scaling for 
h with metallicity holds 



where Z is the metallicity of the gas relative to the Milky Way (i.e., Z = 1 for Milky way 
objects) and ho thus is the width of the C° region in clumps oi Z = 1. We can justify 
this scaling law by the following argument: the abundance of C° is determined by the 
balance between production of C° via photodissociation of CO, and destruction of C° via 
photoionization into C+. The rates for both photoprocesses are dominated by dust shielding 
of the UV photons, while other processes that do not scale with the metallicity (such as 
CO cross-shielding with H2) play only a secondary role. Observationally, we know that the 
dust-to-gas ratio scales approximately as the metallicity (e.g., Koornneef 1982; Bouchet 
et al. 1985), which is not surprising since heavy elements are needed to create grains. 
Therefore the extent of the region over which both photoprocesses balance will be inversely 
proportional to the metallicity. In other words, the extent of the C° region is constant in 
visual extinction (A^), but since the length scale corresponding to A^ scales as 1/Z so does 
h. Note that the CO self-shielding, which depends on the CO abundance and ultimately on 
the C content of the gas, scales in the same manner as the dust shielding. 

Similar considerations lead us to write the scaling law for the extent of the C"*" region 
(/i+) as 



where a is an assumed constant ratio between and h. From chemical modelling (e.g. 
Hollenbach, Takahashi, & Tielens 1991) we surmise that a ~ 1 — 2. 

All the material that is not photoionized (C^) or photodissociated (C°) is in molecular 
form (CO), and given the radius i? of a clump it is easy to compute the size of the CO 
region as /ico = R — h — h^. Figure |l|a illustrates how these scaling laws apply to a sequence 
of clumps with constant radius and decreasing metallicity. Conversely, Figure illustrates 
the scaling laws applied to a sequence of clumps with constant metallicity and decreasing 
radius. 
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From the geometry and scaling laws the contribution of a clump of radius R to the 
number of particles J\f of each species is 



xr(r+\ - l^nijciR-^ - {R - ahy) R > h+ 

■^^^ ^ - 1 IttuvcR' 0< R <h^ 

InnriciiR - ahf - {R - {a + l)hY) R > h + h+ 

Ar(C°) = <; lTTnr]clR - ahf h+ < R < h + h+ (4) 

0< R <h+ 

\r(nn\ _ j t^nric{R - {a + l)h)^ R > h + h+ 

^^^^^ - 0< R <h + h^ 



where rjc is the carbon abundance and we are implicitly assuming that all the carbon 
is in the form of the dominant species in each zone. The abundance of carbon scales 
proportionally to Z, that is 



r/c = VoZ (6) 

where ?7o is the carbon abundance at a metallicity Z = 1. 

We will now compute the mean contribution of a clump to the column density of each 
species, using CO as an example. To eliminate the dependence on clump size, we integrate 
Equation |^ over the clump size distribution R~^ to obtain the mean contribution of a clump 
to the number of CO molecules in the telescope beam. 



1 f J^max 

Ar(CO) = — / M{CO)R-^dR (7) 

where Rmin and Rmax are the lower and upper limits of the clump size distribution, and A4 
is the normalization integral 

R-^dR (8) 

To simplify the calculations, we relate the lower and upper end of the clump size distribution 
to /lo by defining Rmin = 7^o and Rm,ax = ^^o, where 7 and 6 are unitless. 

The total number of particles in the beam will be the number of clumps in it times the 
expected contribution per clump J\f (CO). Using /5 = 3 and the fact that 5 ^ 7 we obtain 
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an algebraic expression for A/"(CO) which we convert to column density (and somewhat 
simplify) by recognizing that the expectation value for the radius of a clump is 



R = -"^^^krT^-TTTT = 27/^0 (9) 



Thus, dividing A/'(CO) by tcR we obtain the mean CO column density per clump 



No 







(a+l)3-6(a+l)^<5Z+3(a+l)^^Z^+253z3 
2(52^2 

+3(a + l)log(^) 



Z > 



g+l 



(10) 



where we define A^'o = |^^H^o^o• Notice that when Rmax < h + h+ even the biggest clump 
is completely photodissociated. As a result there is a critical metallicity Zqo = {I + a)/S 
below which CO is absent from the ensemble of clumps. 

Proceeding similarly for Eqs. § and § we obtain the following expressions for the mean 
C"*" and C*^ column densities 



Nr 



No 
No 5Z 

No 



252^2 



+ 3alog(f ) 



No 




(l+3a+3a2)+6(l+2a)(5Z-3(?2z2 
252^2 

+3«log(^)+31og(^ 

25^25 h -3a log y-^J 



f > 



a+1 



> 



z 

z >0 



>f 



z > 
z > 



g+l 
5 



j> z >0 



(11) 



(12) 



Equations [Tl| and |T2] constitute model A. 



Note that in this parametrization a, the ratio of the sizes of the C~^ and C° regions, is 
also approximately the ratio of the C^ to C° columns times a geometrical correction factor 
of order unity. This stems from the construction of the model (a is the ratio between the 
extents of the C"*" and C'' regions) as well as from the limit of Equations |Tl] and |T2| when 
— ^ 1 and a ~ 1 (c.f., [P). In this limit the ratio Nc+ / Nco is approximately 



Nc+ ^ 31og(f )-1.5 

Nco 



31og(f^ 



4.3 



a 



(13) 
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The value of S determines mainly the CO column, since the biggest clumps are mostly CO. 
Consequently a decrease in 6 will increase both the N^^^+fNco and the N^^^o/Nco ratios. 
Thus the effects of both parameters are relatively well separated. 



3.2. Model B 

We will reconsider now the scaling law for the C° region, expressed in Equation |1]. If 
most of the C° is not originated by photoprocesses but by other chemical reactions not 
strongly dependent on UV photons, the neutral carbon column density may not scale as 

(c.f.. Equation |I|) but as a slower function of Z. As the limiting case, we will assume 
that the size of the C" region is constant with metallicity 



h = h 







(14) 



The size of the C"*" region nevertheless will continue to scale with Z ^, as C"^ is indeed 
produced mainly by photoionization. Therefore Equation still holds. The parameter a is 



then akin to the ratio of C^ to C° column densities at the relative metallicity Z 
corresponding equations for the mean column densities of C+, C" and CO are 



1. The 



Nc+ 



CO 



CO 



No 
No 5Z 

No 



*l?^ + 3«log(f) 



No 


No 




3a^+3a{l-4S)Z+{l-6S+35^)Z^ 

+3«logfe)+3^1og(^) 

-6a^5Z+3a5^Z^+2S^Z^ -U ^n- lr,o- ° ^ 
252^2 '^^'^^^^{SZJ 



Z >f 
f > Z > 



f > Z > 



a3+3a2(l-2i5)Z+3a(l-4^+<52)z2+(l-6i5+3^2+2^3)z3 
252^2 



+3(a + Z)log(^) 



<5-l 



> z >0 



(15) 



(16) 



(17) 



Equations |T5|, |T6| and constitute model B. We compare the relative merits of both 
models in §|^. Notice that models A and B produce identical results at Z = 1. 



4. Computing integrated line ratios 



We can now compute the line intensities using our calculated column densities. Notice 
that the column density ratios between species depend only on two non-dimensional 
parameters: a (the ratio of sizes between the C"*" and regions) and 5 (the ratio between 
the sizes of the biggest clump and the C*^ region). To simplify the calculations we will 
assume optically thin emission in local thermodynamic equilibrium (LTE) for and C°. 
The line integrated intensity (erg s~^ cm~^ sr~^) is then 



hv 



1 n p kTe- 

1= huAui Z^ ' N (18) 

where h and k are Planck's and Boltzmann's constants; u, A^i and T^x are respectively the 
transition frequency, Einstein A coefficient and excitation temperature; Qu is the degeneracy 
of the upper level; Q is the partition function, and is the total column density of the 
species considered. Typically C+ and C° have opacities r < 1 (Zmuidzinas et al. 1988). 
The assumption of LTE will be valid as long as most of the emitting gas is at a density 
greater than the critical density of the transition (~ 2500 cm~^ for [C II], ~ 500 cm^^ for 
the 609 /im transition of [C I]). 

The assumption of optically thin emission from a population in LTE is, however, 
unrealistic for CO. For computing CO (J=l 0) intensities we will use the proportionality 
factor between hydrogen column density and CO line intensity observed at Galactic 
metaUicity, Xco = 3 x 10^° cm-V(K km s'^) (Scoville & Sanders 1987). The CO and 
hydrogen column densities at Z = 1 are related via the CO abundance rjo, Nco = N^riQ. 
Thus, the integrated line intensity (in erg s~^ cm~^ sr~^) can be related to the CO column 
density in the following way 



Ico = (19) 

where Aio is the wavelength of the J=l CO transition, and the factor ^ comes from 
the conversion between velocity-integrated antenna temperature and frequency-integrated 
specific intensity. In this expression we are implicitly assuming that the emissivity of CO, 
determined by the excitation conditions of the molecule, does not depend on Z. All the 
dependence of Ico on the metaUicity stems from the metaUicity dependence of A'co- This 
assumption will be discussed further in §||. 

For the optically thin emission from [C II] (^P3/2 158 /im) and [C I] (^Pi ^^Po, 



609 //m) we assume excitation temperatures of 200 K and 30 K respectively (Wolfire, 
HoUenbach, & Tielens 1989; Wright et al. 1991). Thus, 



I[cii] = 1.28X 10-'^iVc+ 

I[ci] = 7.54 X IQ-^^iVc^ (20) 
Ico = 1.73 X 10-26 ]Vc^ 

where we arc using the canonical value r)o = 3 x 10~^ (e.g., Tielens & HoUenbach 1985a) 
and the line intensities are in units of erg s"^ cm-^ sr"^. The equations for line ratios then 
are 



I[cii]/Ico = 7.4x IC^iVcT/iVTO 

I[ci]/Ico = 4.4x lO^iVc^/iVc^ (21) 
I[cii]/I[ci] = 1.7x lO^iVc^/iVc^ 

where only two of the three equations are independent. 

Note that because [C II] and [C I] are a two and a three level system respectively, the 
intensities are independent of the precise value of the temperature as long as it exceeds the 
line excitation threshold. 



5. The values of the pEirameters a and S 

We can use Galactic data to assign values to a and 5. In particular, we could use 
COBE's FIRAS measurements averaged over the entire galactic plane. Unfortunately 
FIRAS did not obtain a significant detection of the CO (J=l 0) transition. This opens a 
range of possibilities for the CO intensity: 1) the formal value obtained by FIRAS (0.6 ±0.7 
erg s~^ cm~2 sr~^; Wright et al. 1991), 2) the result of a multi-line excitation fit performed 
with the higher transitions of CO which FIRAS did detect (0.5 erg s^^ cm^^ sr~^; Wright 
et al. 1991), or 3) the result of assuming the standard empirical ratio I[cii]/Ico — 4400 
(Wolfire et al. 1989), which corresponds to a CO (J=l — > 0) intensity of 0.4 erg s~^ cm"^ 
sr~^ given the observed FIRAS [C II] intensity of 1726 ± 70 erg s~^ cm^^ sr~^. Adopting 
the intermediate value Ico = 0.5 erg s~^ cm"^ sr~^ we obtain the ratios I[cii]/I[ci] — 270 and 
I[ci]/Ico — 13 from the FIRAS measurements. Numerically solving the system of Equations 
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21 we find the values ol = 1.25 and 5 = 435. Note that this result is independent of the 



choice of model A or B, since as was pointed out before both models are identical at Z = 1. 

The resulting column density and line intensity ratios are shown in Figure ^. With 
the above values for a and 6, the ratio of C+ to C'^ column densities in model A yields 
Nc+/Nco ^ 1.65 ± 0.10 over a wide range of metallicities. In model B the quantity that is 
approximately conserved for different metallicities is the ratio of C° to CO column densities, 
and consequently the ratio of intensities. 

It is interesting to point out again that lowering the value of the maximum clump 
size Rmax (set by 6) increases the column densities of C"*" and C°. This occurs because 
the biggest clumps contribute mostly molecular (CO) gas. In fact, they dominate the CO 
column density. The effects on the line ratios of varying the maximum size parameter 6 
are shown in Figure ^ Some of the abnormally high [C II] /CO intensity ratios observed in 
certain sources could be due to a local change in the clump size distribution, caused by the 
disruption of the biggest clumps. 



6. The effects of density and radiation fields 

Up to this point, our parametrization has been purely geometrical. We have not 
explicitly included changes in the gas density (n) or the radiation field intensity (xuv)- We 
can introduce these changes through their effects on the parameters a, 6 and Hq. 

According to chemistry calculations (van Dishoeck & Black 1988), at a hydrogen 
density n = 10^ cm~^ the C"*" column density increases by a factor of ~ 10 when Xnv is 
increased from 1 to 10^ (in units of 1.6 x 10~^ erg cm~^ s~^; Habing 1968), while the C° 
column density increases roughly a factor of ~ 2 to 3. Thus, the column density of C° 
(controlled by the parameter Hq) as well as the ratio of the C^ to C° column densities 
(controlled by the parameter a) are both somewhat affected by the UV field strength, albeit 
not dramatically. The logarithmic dependence of the column densities on the radiation 
field occurs because, at densities larger than ~ 10^ cm~^, the column density of the species 
produced by photoprocesses varies proportionally to the penetration depth of the UV 
radiation, which in turn increases logarithmically with Xuv (Stacey et al. 1991; see also 
discussion in Mochizuki et al. 1996). On these grounds, we surmise that both a and ho are 
directly proportional to the logarithm of Xuv Assuming that the upper end of the clump 
size distribution, Rmax = Sho, remains unaffected by the UV field strength, 5 is inversely 
proportional to Hq. 

A qualitative argument can be used to understand how changes in the density will 
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affect ho, a and S. Higlier densities will yield higher recombination rates for C"*" and the 
more active chemistry will produce more CO at the expense of C*^ and C"*". Although it is 
a priori unclear how the column densities of C+ and C° depend on the density, chemical 
modelling (Ingalls 1998) shows that it is also approximately logarithmic. 

Figure ^ shows the results of the steady-state chemical network calculations on an 
Ay > 8 cloud. Superimposed are the 2-dimensional fits to the chemistry results: 



logNco = 0. 13 log Xuv - 0.18 log n + 18.32 (22) 

logf^^) ^ 0.051ogXuv + 0.021ogn-0.16 (23) 
V Nco J 

These equations are valid over a wide range of physical conditions (xuv ~ 50 — lO''; 
n ~ 10^ — 10^ cm~^) and the errors introduced by using the fits over the actual chemical 
network results are ~ 10%. 



The column densities of and C° obtained from Eqs. |2^ and p3| can be directly 



related to the parameters a and Hq by observing that, to first order, Nqo ^ nrj^hQ and 
Nq+/Nqo ~ a. These two equations, combined with 5 oc hQ~^ lead to: 



h0.05 



a ^ 1.25(^1 (24) 



435 "'^V^V'-'' (25) 



-0.13 

where Xuvo and are the average radiation field and volume density of the regions that 
dominate the FIRAS line emission used to assign a = 1.25 and 6 = 435. Note that because 
HII regions probably dominate the FIRAS emission (Wright et al. 1991; Bennett et al. 
1994), Xuvo ^ 1- Since the FIRAS ratios are similar to those observed in many Galactic 
star-forming regions (in particular Orion) we normalize all values of density and UV field 
to those derived for Orion, namely no ~ 10^ cm~^ and Xuvo ~ 250 (Tielens & Hollenbach 
1985b; Stacey et al. 1993). 



7. Comparison with observations 



The predictions of models A and B can be directly compared with observational data. 
To properly accomplish this comparison we need [C II] and [C I] observations on a sample 
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of objects of the same class with different metalhcities. Unfortunately very little [C I] data 
is available on extragalactic sources of differing metalhcities (c.f., Wilson 1997; Table 2) and 
only a few of them have been observed in [C II]. 

The data we use for the comparison are listed in Table |l|. The selected sources are a 
set of star-forming regions with different metalhcities which have been observed in [C II] 
and [C I]. To them we have added N 27 in the Small Magellanic Cloud (SMC), a source 
that does not yet have [C I] observations but which we expect will be important to probe 
our models because it has the lowest metallicity in the sample. 

Figure ^ compares the I[cii]/Ico intensity ratio predictions of models A and B with 
the observational data, using the values for a and 5 derived from FIRAS measurements. 
Both models give essentially identical predictions for the decreasing trend in the I[cii]/Ico 
intensity ratio with increasing metallicity. This similarity between models A and B should 
not be surprising, since the effect of keeping h constant in a clump (as model B does) is to 
enlarge its CO region by only a small fraction for Z < 1, with respect to the same clump 
in model A. The fact that the I[cii]/Ico ratio for most sources (except 30 Doradus) falls 
along the model line can be interpreted as showing that N 27, IC 10-SE, N 159 and Orion 
all share similar values of n and Xuv Another possibility is that the density and UV field 
in these regions conspire in such a way as to keep their line ratios close to the model line. 
This can happen because, as shown in Equation the effects of density and radiation 
field go in opposite directions and have the potential to cancel each other. Denser regions, 
which are expected to have lower column densities of C"*" and C°, will be in general more 
efficient at producing stars and thus will have on average a higher UV field which in turn 
will increase Nc+ and Nqo . However, because of the large difference in the exponents in 
Equation this cancellation is unlikely to occur and we can safely assume that these 
star-forming regions posses similar densities and radiation fields. 

The large I[cii]/Ico ratio observed in 30 Doradus suggests either that the gas in this 
source is under quite different physical conditions, or that Rmax, the upper hmit of the 
clump size distribution, is actually smaller than for the rest of the sample. Lowering the 
size of the biggest clump (by making 6 smaller) will produce more photoionized material 
at the expense of CO, which will in turn raise the observed line ratio (c.f., §|^). The gas 
in 30 Doradus is under an extremely intense (and hard) UV field, Xnv — 3500 (Poglitsch 
et al. 1995). From the numerical exponents of the relationships postulated in §^ (c.f.. 
Equations ^ and an increase by a factor of 15 in Xm with respect to Orion (c.f.. Table 
[ip is not enough to explain the observed line ratio. The average gas density must also 
decrease to about 2000 cm~^, about 5 times lower than Orion's. This density is within the 
range allowed by previous determinations (Poghtsch et al. 1995). Thus, it is not necessary 
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to invoke a smaller maximum clump size to explain the observed [C II] /CO line ratio, 
although it remains a logical possibility (e.g., Pak et al. 1998). The 30 Doradus nebula is 
an especially active massive star-forming region, and as such it is the site of several violent 
events which could disrupt the larger molecular clumps. 

Figure H illustrates the comparison between the I[ci]/Ico intensity ratio predictions 
of models A and B and the observations. In this case the predictions of the two models 
differ drastically. While model A predicts a ratio that rapidly increases for decreasing 
metallicities, such a trend is not apparent in the observational data, which tends to cluster 
around a ratio of ~ 10, essentially the prediction of model B. Both Large Magellanic 
Cloud (LMC) points obtained by AST/RO (Stark et al. 1997) are lower limits to the ratio 
because of pointing uncertainties. They are roughly consistent with either model at the 
metallicity of the LMC, although N 159 favors model A while the discrepancy in the case of 
30 Doradus is smaller with model B. A stronger constraint is provided by the IC 10-SE [C I] 
observations (Bolatto et al. 1998), which would have to be a factor of ~ 2 more intense in 
[C I] in order to be compatible with model A. An added complexity in the interpretation 
of the ratios in IC 10-SE is that the [C II] observations have a much larger beam size than 
the [C I] (approximately 60 arcsec and 10 arcsec respectively). Therefore, the I[cii]/Ico 
ratio may be sampling parcels of gas under much different physical conditions from those 
sampled by the I[ci]/Ic05 a fact that could make the ratio observed in IC 10-SE appear 
artificially low. 

Stronger constraints on the models would be provided by [C I] observations of the 
SMC. If this trend of constant [C I]/CO intensity ratio is confirmed, we are led to conclude 
that most of the in the clumps is not created by photodissociation but by other chemical 
processes. Otherwise, we would expect the trends of [C II] /CO and [C I] /CO intensity 
ratios with metallicity to be similar, as in model A. 

8. The value of Xco 

The proportionality between the CO emission intensity and the column density of II2, 
A^H2 = -^cqIco! has been observationally established using different methods (see Magnani 
& Onello 1995 for a review of the methods and their caveats). The proportionality is well 
accepted for giant molecular clouds (GMCs), although the precise mechanisms by which 
Xqo becomes relatively independent of the physical conditions in the cloud remain unclear 
(e.g., McKee 1989). 

There are a few observational determinations in the literature of the dependence of Xqq 
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on metallicity Z. Using interferometric CO observations to obtain virial mass estimates of 
extragalactic GMC complexes, Wilson (1995) found a power law dependence with a slope 
of —0.67. Using a combination of single-dish and interferometric CO data and virial mass 
estimates, Arimoto et al. (1996) find a slope of —0.8. 

In the context of the models presented here, we can use Equation |l^ to investigate the 
expected dependence of Xco on metallicity. Accordingly 



^co(^) ^ iVco(Z = 1) 
Xco(^ = l) Arco(^) ^ ^ 

and thus Xqq is inversely proportional to the mean column density of CO. The resulting 
slope with metallicity, according to Equations |1^ or |1^, is ~ —1. 

A comparison of our models with the results from Wilson (1995) and Arimoto et al. 
(1996) is shown in Fig. |^. The precise value of the clump size distribution exponent /5 or 
the parameters a and b do not affect the value of this slope, which mainly depends on the 
dominant term of the polynomial part of Nqo- Although the Wilson (1995) data points 
certainly fit a —0.67 slope better, the quality of the data does not exclude a slope closer to 
— 1. A similar situation occurs with the Arimoto et al. (1996) data. 

One of the reasons why the actual slope could be different from the theoretical 
expectation is a systematic change in the CO excitation temperature with metallicity. Such 
a dependence might be expected, since both the heating (F) and cooling (A) of the gas are 
affected by metallicity. The F term is fundamentally controlled by photoelectric ejection of 
electrons from dust grains, which collisionally heat the gas. The cooling term is controlled 
by emission in the [C II], [O I] and [C I] fine structure transitions. The balance of the 
system at different metallicities and dust-to-gas ratios has been computed by Wolfire et al. 
(1995), who conclude that the temperature of the gas is higher for low Z . If the excitation 
temperature of CO is higher at low metallicities, the column of CO needed to produce 
a given intensity Ico will drop, together with the corresponding column density of H2. 
Consequently, there will be a drop in the Xco factor as Xco = A^H2/Ico which will tend to 
fiatten the curve of Xqq,{Z\ 



9. Summary and conclusions 

We have presented a model of PDRs of varying metallicity that depends on only two 
free parameters (a, the ratio of sizes between the C^ and C'' regions and 5, a measure 
of the size of the biggest clump in the ensemble). This approach naturally incorporates 
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the dumpiness of the interstellar medium into the calculations and, when combined with 

some simple scaling laws, allows us to study the influence of metallicity on the PDR line 
emission. We have shown how to modify it to include the effects of density and radiation 
field, and we have used it to study the dependence of the l[cii]/Ico and l[ci]/Ico line ratios 
on the metallicity. In comparing our calculations with observational data, we find that the 
model explains well the observed trend of enhanced [C II]/CO line ratio with decreasing 
metaUicity as the natural result of the augmented fraction of photodissociated gas in a 
clump. We have shown that enhanced ratios can be produced by lowering the density 
or increasing the radiation field in the ISM, but also by lowering the upper limit of the 
clump size distribution as may happen in 30 Doradus (Pak et al. 1998). We find that 
the available data favors a [C I] /CO intensity ratio essentially independent of metallicity, 
albeit the paucity of observations does not exclude other possibilities. This is difficult to 
understand if most of the C° is produced by UV photons as is the case with C+, arguing for 
a chemical origin for at least part of the neutral carbon in these star-forming regions. More 
[C 11] and [C I] observations on low metalUcity sources such as the Small Magellanic Cloud 
are necessary, however, to establish this point. Finally, we have studied the prediction of 
the model for the trend of the Xqo factor with metallicity and compared it to previous 
observational studies, and we find that the agreement is reasonable. 
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Table 1. Source metallicities and intensity ratios 



Source 


12+log(0/H) 


Xuv 


n 


IlCiil/Ico" 


I[ci]/Ico" 


References 








(cm-3) 




Milky Way 


(8.75) 






3450 " 


13 " 


Wright ct al. 1991, Bennet et al. 1994 


Orion 


8.75 


250 


10" 


4500 


10 - 15 


Lcqucux ct al. 1979, Stacey et al. 1993, 
Tauber et al. 1995 


N 159 (LMC) 


8.43 


300 


8 X 10^ 


5600 


> 20 


Dufour 1984, Israel et al. 1996, 
Stark et al. 1997 


30 Doradus (LMC) 


8.43 


3500 


103 - 10* 


69000 


> 8 


Poglitsch et al. 1995, Stark et al. 1997 


IC 10 SE (IC 10) 


8.17 






14000 


19 " 


Lcqucux ct al. 1979, 

Madden ct al. 1997, Bolatto et al. 1998 


N 27 (SMC) 


8.02 






13000 




Dufour 1984, Israel et al. 1993 



Rubio et al. 1993 



^Intensities in erg cm ^ s ^ . 
''Assumed to be identical to Orion's. 
''Computed using Ico=0-5 erg cm~^ s~^. 

''Also from multiline excitation study using Johansson et al. (1994) observations. 
« Error- weighted average for points SL2, MCI, MC2 and MC3. 



-22- 




Fig. 1. — (a) Top Row. A sequence of clumps with identical radius and decreasing metallicity 
according to model A. When Z — {a-\- l)ho/R, the CO is completely photodissociated and 
the clump consists only of C° and C+. The next stage happens at Z — aho/R, when all the 
neutral carbon is photoionized into C"*". At even lower metallicities all the gas consists only 
of C"*". (b) Bottom Row. A sequence of clumps with constant metallicity and decreasing radii 
in model A. Similar sequences where the limiting cases occur at slightly different metallicities 
can be created using model B. 
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Relative Metallicity (Z) 

Fig. 2. — (a) Top. The mean column densities of CO, and C"*" versus metallicity resulting 
from models A and B. The values used for the hydrogen density and the size of the region 
for the purpose of the plot are n — 10^ cm~^ and ho — 10^^ cm (~ 0.03 pc). (b) Bottom. 
The line intensity ratios versus metaUicity for models A and B. Notice that, in model A, 
the [C 11] /[C I] ratio is almost constant over a wide range of relative metallicities, while in 
model B the [C I] to CO ratio is approximately constant. 
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10' 10^ 10^ 10^ 

Largest clump size (6) Largest clump size (5) 

Fig. 3. — The influence of the upper hmit of the clump size distribution (i.e., the parameter 
S) on the model ratios at Galactic metaUicity {Z — 1). The biggest clumps contribute mostly 
CO to the ensemble; consequently the line ratios decrease for higher clump sizes. 
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2 3 4 2 3 4 



log iXyjy) log (z^^) 

Fig. 4. — (a) Left Column. Two dimensional fits to the results of a chemical model (Ingalls 
1998) for logNcjo and \og{Nc+ /Nco). The gray scale corresponds to the results of the 
chemical network. The black contours are traced at the same levels and show the quality of 
the fits. (Top Left) Contour levels are IQ-*^^ to IQ-*^^ cm~^ with logarithmic spacing 0.2. 
(Bottom Left) Contour levels are -0.05 to 0.1 with logarithmic spacing 0.05. (b) Right 
Column. Errors introduced by using the fits (F) rather than the actual results of the chemical 
network (A'"), expressed as \N — F\/N. The RMS error for Nqo and Nq+INqo over the whole 
range of density and UV field shown here are 11% and 8% respectively. Contour levels are 
10% to 30% by 10%. 
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7.6 7.8 8 8.2 8.4 8.6 8.8 

12 + log(0/H) 



Fig. 5. — Comparison of model A and B predictions with [C II] /CO intensity ratios observed 
in several sources. The solid lines are the calculated [C II] /CO (J=l —>■ 0) intensity ratio 
with a and 6 derived from the COBE FIRAS ratios. The dotted line corresponds to the 
COBE-based a and 6 corrected in the way prescribed in by assuming that for 30 Doradus 
Xnv = 3500 and n = 2000 cm~^. The data points are: the Orion interface region (triangle), 
N 159 in the LMC (circle A), 30 Doradus (circle B), IC 10 SE (square), and N 27 in the 
SMC (star). The error bars in the ratio are calculated assuming 30 % calibration accuracy 
in both lines. 
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Fig. 6. — Same as figure |^ but for [C I]/CO intensity ratios versus metallicity. Tlie point for 
IC 10-SE correspond to tlie average ratio in tfie complex (Bolatto et al. 1998). The point 
for 30 Doradus has been shghtly displaced in metallicity to eliminate confusion in the error 
bars. 
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Fig. 7. — The CO (J=l 0) intensity to hydrogen column density conversion factor (Xco) 
versus metaUicity, as predicted by the models and deduced from observations. The data 
points are those used by Wilson (1995). The solid hne corresponds to Wilson's best fit with 
a slope of —0.67. The dotted hne is a similar result obtained by Arimoto et al. (1996) using 
an independent data set. The slope of Arimoto's best fit is —0.8. The dashed line represents 
the result of the models here developed, as explained elsewhere in the text. The slope is 
—1.05. Models A and B produce nearly identical results. 



